Coherence and clock shifts in ultracold Fermi gases with resonant interactions 
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Using arguments based on sum rules, we derive a general result for the average shifts of rf lines in 
Fermi gases in terms of interatomic interaction strengths and two-particle correlation functions. We 
show that, near an interaction resonance, shifts vary inversely with the atomic scattering length, 
rather than linearly as in dilute gases, thus accounting for the experimental observation that clock 
shifts remain finite at Feshbach resonances. 
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Interatomic interactions limit the accuracy of atomic 
clocks, causing density-dependent clock shifts in radio 
frequency (rf) transitions. Similarly, such shifts play an 
important role in probing correlations in atomic gases [H , 
where, e.g., rf spectroscopy has been used to detect the 
presence of molecules and provide evidence for pairing 
gaps did, 3- Surprisingly, experimentally observed clock 
shifts become small when interactions are resonantly en- 
hanced [1,0], a result we explain here, first by develop- 
ing a general theory of the average clock shift and then 
showing that in the strongly interacting regime the shifts 
depend inversely on interatomic scattering lengths. 

Single component spin-polarized Fermi gases do not 
experience clock shifts, since the rf couplingpreserves 
polarization, forbidding s-wave interactions [5[. In ad- 
dition, in mixtures of interacting fermions in two states 
|1) and 1 2), e.g., the lowest two hyperfine states of ^Li, 
clock shifts are absent for transitions |1) |2), since the 
interaction energy is invariant under the rf field In- 
teractions in two-state mixtures can be probed, rather, 
by driving transitions to an initially empty state, e.g., in 
^Li, from |2) to the next hyperfine state |3). 

Rf transitions in atoms are usually described in terms 
of coherent evolution of a two-level system undergoing 
Rabi oscillations, represented by rotations of an equiv- 
alent pseudospin on the Bloch sphere. In the present 
problem, the rf field rotates an atom in state |2) into a co- 
herent superposition, |2) |/?) = cos0|2) -I- e"*"* sin6'|3). 
(The angles 4> and 9 depend on time, the rf pulse power, 
and its detuning from resonance.) To the extent that in- 
teratomic interactions only shift the energy levels, but do 
not broaden the lines, the Bloch sphere picture is valid. 
However, strong interactions lead to an incoherent and 
irreversible evolution, and in the long-time, weak pulse 
regime the probability to find the system in a particular 
final state is given rather by Fermi's Golden Rule. 

We base our discussion on linear response theory (in 
Ref. Q fewer than 30% of the atoms are transferred 



by the rf pulse, with comparable numbers in Ref. 0]). 
We consider a spatially uniform Fermi gas with three 
internal states, and Hamiltonian H ~ Hq + Hs, where 
Hq ~ X]?=i ^i-^ii £i is the hyperfine plus Zeeman energy 
of level and Ni is the total number of atoms in state 
The system Hamiltonian is (?i = 1 throughout), 



H, = J d\ (^i-VV4(r) • VV',(r)' 
+ I d\d\'v.,{v-v')4{v)i,]{v')^l;,{v')i,.,{v). (1) 



Here Wij(r — r') is the bare potential between two atoms 
in \i) and |j) (not a low-energy effective interaction). 

The (essentially spatially uniform) rf field primarily 
couples the states |2) and |3); other internal transitions 
are far off-resonance. We therefore describe the rf cou- 
pling by iJit — B{t)Y with the pscudospin-flip operator 

Y^i j d\{i;l (r) V^2 (r) - 4 {v)^^ (r)) , (2) 

where B{t) ~ 2lj^ cos ujI, with the Rabi frequency. 

We assume an initial many-body state |12) of energy 
Ei2, containing atoms in states |1) and |2). The transi- 
tion rate of atoms from state |2) to |3), I{oj), is 

I{uj) = 27r^i^ |(/|y|12)p^(c.+£;i2-£;/) = 2ujix"{^), 

f 

(3) 

where the sum runs over all eigenstates |/). The mean 
frequency for the transition. 



where = £3 ~ £2, and Qc is the average clock shift. Us- 
ing Eq. ([3]) we derive the f-sum rule for the first moment 
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of the spectrum, 

/OO 7 
—cox"{co) - Y.'^Ef - Er2)\{f\Y\12)\ 
-OO j: 



-(i2|[[y,i/],y]|i2). (5) 



For state |3) mitially unoccupied, 

/ — X {^) 
dVdV'(12|v4(r)^3(r)V^^(r')V'2(r')|12) -iV2, (6) 



where N2 is the number of atoms in state |2). We thus 
derive the simple result for the mean clock shift 



^^c = ^(i2|[[r,i?J,r]|i2) 



(7) 



More generally, the expectation value becomes an ensem- 
ble average. 

The average shift is simply the change in energy per 
2-atom required to rotate all |2) atoms in state |12) in- 
finitesimally into state |/?). This rotation transforms 
1 12) into a state = e^''^|12) of |1) and /3 atoms, 

with the same spatial many-particle wave function as 
1 12). The energy difference of the two states is thus 
5E = (l/3|i?|l/3) - (12|H|12) ^ ie'2(12|[[r,i7],y]|12) 
to second order in 9. Since the number of atoms 5N 
transferred from |2) to |3) is N2&^ ■, the change in system 
energy per atom is given by Eq. ([7]), a result valid both 
for superfluid and normal phases. The general expression 
for the clock shift in terms of the interactions is 

n,^ — f d3r[i;i3(r)~z;i2(r)](Vl(r)V|(0)V2(0)Vi(r)), 
"2 J 

(8) 

where n2 = N2/V, V is the system volume, and (...) 
denotes an ensemble average. 

To calculate the clock shift explicitly, we assume for 
simplicity in the following that the interatomic potentials 
arc short range contact interactions, each with a large 
momentum cutoff. We also assume that ni = n2 = n. 
The bare couplings, gij, are related to the measured low 
energy couplings gtj = 47ra,;j/m, where atj is the scat- 
tering length, by g:r.^ = g^^ + d^p/{2Tr)^{m/p'^) = 
g~^ + mAij/27r^, as one sees by solving the usual t- 
matrix equation. Thus^y = 47r/ (m(a~-^ -I- r'^l^f), where 
fo.ij = 7r/2Aij. Physically we expect rg.ij to be of order 
the characteristic length [CQm/m^Y/^aQ for the inter- 
atomic potential, where Cg is the strength of the van der 
Waals interaction; in ®Li, rg ~ 6809. For short range 
contact interactions, Eq. ([S]) becomes 



-(513 - ffi2)(^I(0)^|(0)^2(0)Vi(0)). (9) 



The correlation function is related to the free energy, 
F , of the system by the thermodynamic identity 

^^EM. ^ ^{dH/d-g,2) - (V^I(0)V;t (0)^2(0)^1(0)). 

(10) 

We now take tq to be independent of the states of the 
atoms. From Eq. ^ we then find the mean shift: 



(-) 






\al2J 


V.913 


512/ 



1 d{F/V) 



[ID 



^12 



When the magnitude of the scattering lengths are 
^ ro, which for ^Li is always valid above B = 600G, 
the factor 013/(112 is ~ 1; the shift then involves only 
the renormalized interactions, depending inversely on the 
scattering lengths. This result has the same qualitative 
behavior, in the regime when k\a\ 3> 1, where fc is a 
typical particle momentum, as predictions based on a 
mean-field shift calculated from the real part of the for- 
ward scattering amplitude in a dilute gas in the absence 
of effective range contributions Q . 

The derivative d{F/V)/dgi2 < behaves perfectly 
smoothly at the resonance in the 12-channel, and at zero 
temperature its value there is ^ —Cp' where ep = p|/2m 
is the free particle Fermi energy and pf is the Fermi mo- 
mentum. (At characteristic densities ^ 2 x lO^^cm^'^, 
Pp^ ~ 2000ao |3|.) We assume that ni = 712 = n 
from here on. From the Monte Carlo calculations in 
Ref. of the energy E at zero temperature we estimate 
that d{E/V)/dgi2 — — 0.50e| at the 12-resonancc at 

B ~ 834 G. Then = giid{F/V)/dg^2 Pp/mais; 

this result, valid for 013 ^ tq, scales as ^/ep. On 
the other hand, at a 13- resonance {B ~ 690 G), fie = 
— l/{ngi2)d[F/V)/dg^2- s^^j the average clock 

shift remains finite, even when scattering is resonant. 

Figure 1 shows the correlation function d{E/V)/da^2 
calculated from the data of Ref. Q as a function of B 
for equal populations at zero temperature. In the BCS 
limit, the correlation function tends to the Hartree-BCS 
result, + A^/5^21 where A is the gap, while in the 
BEG limit (pFai2 ^ 1) the system is composed en- 
tirely of 12-molccules with binding energy l/mai2, and 
dE/dgi2 — — 47™/m^ai2, diverging as ai2 0. For 
small positive scattering length in a gaseous atomic state 
limit the correlation function tends to the Hartree value 
72^, as on the BGS side. The large drop in dE/dai2 seen 
in Fig. 1 for B below 700 G reflects the presence of 
12-moleculcs in the initial state. 

In the limit in which the bare couplings are small, |a| <C 
To, the bare and renormalized couplings are equal, and 
the correlation function becomes simply ; thus 



^^c = ?^(.gi3 - 512), 



(12) 



712 



the expected weak coupling result. The intermediate 
regime in which tq <^ \a\ <C p^^ is more complicated 
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FIG. 1: Correlation function d{E/N)/da-^ in units of p-p /2m 
calculated from the data of Ref. 



since here short range correlations as well as bound states 
in the appropriate channel play a role. From Eq. (jlip 
we find in this regime for an initial state without 12- 
moleculcs that 
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(13) 



which differs from the weak coupling result p2|) by a fac- 
tor 912/913- To derive this result, we note that the cor- 
relation function in the absence of molecules is given by 
the Hartree approximation, n^, multiplied by the square 
of the Jastrow factor that takes into account short-range 
correlations due to the interatomic interaction. At low 
energy, the Jastrow factor at the range of the forces is 
essentially 1 — 012/ tq ~ .912/512, the ratio of the exact 
to the non-interacting two particle wave function. Thus 

{^plio)4io)Mo)Mo)) = (512/512)^^2 [HI. 

The difference between the result (fT3)) and Eq. p2)) lies 
in the distribution of spectral weight for transitions in the 
two cases. The first term gi^ — gi2 in the last expression 
in Eq. is the weak coupling result. 

The second term however arises from short range cor- 
relations, including, for (713 > 0, 13-molecules. The 
squared matrix element for free-bound transitions is pro- 
portional to {913 — 912)'^ 0, and the l/(?i3 dependence 
reflects the fact that the binding energy of the molecule 
{l/maf^) increases with decreasing 313. (Strongly bound 
13- molecular states are not described by the simple con- 
tact interaction model, but frequencies of transitions to 
them in general lie outside the range explored in exper- 
iment.) For ai3 < 0, there are no bound states, but 
higher- lying states with kinetic energies of order l/ma^3, 
arising from short-range correlations, extend the spec- 
tral weight out to frequencies of order ~ l/mai^, and 
contribute ^ — l/ai3 to the sum rule [l^ . 

In the BEC limit, pFai2 <C 1 the system is com- 
posed entirely of molecules of binding energy Eb ~ 
1/7710^2 ■ The energy density of the system is just E/V ~ 



i/ma^2j and from Eq. (|lip we obtain 



Qc = 2Eb{1 - ai2/ai3), 



(14) 



as can also be calculated from Ref. [91 . 

In applying our results to experiment it is important 
to consider the extent to which, for ai2 > 0, molecules 
are present in the initial state. The experiment of Ref. Q 
is carried out on a timescale short compared with that to 
create diatomic molecules via three-body collisions. Thus 
the correlation function entering the sum rule is that of 
the molecule-free metastable state. On the other hand, 
in the initial state in the experiment of Ref. Q molec- 
ular states were equilibrated, and their effects must be 
included in the correlation function. Furthermore com- 
parison with experiment requires that measurements be 
carried out over the entire frequency range, on scales 
~ 1/7710^3, where there are significant contributions to 
the spectral weight. 

The solid line (a) and (b) in Fig. 2 shows the predicted 
clock shift for equilibrated ^Li, in units of ep/h in the 
region ~ 600-1200 G at zero temperature. Scattering 
lengths are taken from Ref. []~3|, and we assume pp ~ 
l/2000ao \ The dashed-dot line (c) shows the result (fM)) 
for a non-interacting gas of molecules. The dashed part of 
curve (a) below the 12-resonance, the difference of the full 
Monte-Carlo result (b) in the 12-molecular region, and 
the two-body contribution (c) , is a measure of the many- 
body contribution to the shift below the 12-resonance. 
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FIG. 2; Calculated clock shift for a gas of Li atoms at zero 
temperature as a function of magnetic field. See text for de- 
tails. 

The correct picture of rf field excitation combines both 
coherent rotation of the initial state with transitions that 
change the number of elementary excitations of the sys- 
tem and thus lead to a width, F, of the spectrum I (to). 
For times short compared to uj^^ and F^^ the rf field 
indeed starts to rotate all |2) atoms coherently into |3). 
For a strong drive, with lob. ^ F, the many-body state 
undergoes Rabi oscillations, damped on a time scale 
^ p_i^ rj,j_^^ long-time behavior, when all oscillations 
have damped out, is captured by the Golden Rule. 
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Under an instantaneous rotation of |2) to the \lf3) 
state produced has the same spatial wave function as 
1 12). However the bare amphtude for 1,(3 —* l,/3 scat- 
tering, gip = cos^ 6'.9i2 + sin^05i3, differs from the bare 
12 amphtude by 6g = gip - gu = sin^ 9{gi3 - §12). 
One can understand Eq. ^ from this point of view, 
by noting that since the spatial wave functions of |l/3) 
and 1 12) are identical, the difference of their energies 
(and free energies at finite temperature) is, for small 6, 

Eip ~ Eu = Sg{^liO)4iO)MO)MO)). 

The rotated state |l/3) is not in fact an energy eigen- 
state, and thus the response has a width. Under rotation 
the interaction energy becomes 

^^int = [gi2 COS^ e + Sin^ e]ll>\ll)'^pi!j3'4)i 

+ [513 cos^ e + 512 sin^ 6']V'IV'q'0q'0i 
+ (513 - 512) sin 6* COS 6* i'\ (ijjli'p + ipl'ipa'j i'l, 

(15) 

where |a) = cos 6*1 3) — sin 6*12) is the superposition of sin- 
gle particle states |2) and |3) orthogonal to |/3). The first 
term in Eq. (jl5[) is the original interaction with a modi- 
fied coupling constant, in terms of states |1) and |/3), in- 
stead of |1) and 1 2). The second term acting on a rotated 
state without a atoms present vanishes. The final term, 
however, mixes the two particles in states |1) and |/3) into 
1 1) and \a) . Thus while the external rf field coherently ro- 
tates the initial state, the latter interaction decoheres the 
rotated state, leading to a mixture of energy eigcnstates 
and a response with nonzero width. Even in the absence 
of the final term, the rotated state is not an eigenstate of 
Hg, since the spatial wave function of an eigenstate with 
gip differs from that with the original gi2. 

The lowest order perturbative approach does not in- 
clude collective effects, e.g., screening of the interaction 
by the medium, and fails to satisfy the sum rule in the 
presence of interactions. As indicated in Ref. Q, these 
effects must be included by summing chains of bubble 
diagrams. Such terms are especially important in the 
superfluid state, since they take into account transitions 
to states in which the phase of the condensate is altered 
globally. 

In summary, we have shown from sum-rule arguments 
that the average clock shifts in strongly interacting Fermi 
gases are finite at Feshbach resonances. The interaction 
model we have employed is admittedly crude, since it 
does not allow for the many bound states that exist for 
real interatomic potentials. Detailed calculations of the 
spectral distribution must be done microscopically. In 
addition to the self-energy corrections to propagators, it 
is necessary to include vertex corrections, as pointed out 
in Ref. Q • Such effects are important even in a low den- 
sity gas, where one finds interference terms between 12- 



scattering and 13-scattering processes [13|]. These pro- 
cesses, which correspond to the Aslamazov-Larkin dia- 
grams for fluctuation-induced effects in superconductors 
[l5|, will be discussed in a future publication. In addition, 
for detailed comparison with experiment it is necessary 
to take into account the inhomogeneity of the particle 
density. 
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